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Abstract 

We investigate a ID quantum system associated with the Ising model in a 
field (the dilute A3 model) by the recently developed quantum transfer matrix 
(QTM) approach. A closed set of functional relations is found among variants of 
fusion QTMs which are characterized by skew Young tableaux. These relations are 
proved by using a quantum analogue of Jacobi-Trudi formula, together with special 
features at "root of unity" . The numerical analysis on their eigenvalues shows a 
remarkable coincidence with exponents characteristic to £"8. From these findings, 
we have successfully recovered the Thermodynamic Bethe ansatz equation by 
Bazhanov et al, however, without specific choice of strings solutions. 
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1 Introduction 



Perhaps the most striking prediction from the deformation of conformal field theory to 
statistical mechanics is the hidden Eg universality structure in the Ising model in a field 

0]. This triggers active researches in both field theory]^]- and statistical mechanics 
lMl"[IIH- ^^^^ report, we shall attack the problem in view of a solvable lattice model . 

In a L— state RSOS model is proposed of which Boltzmann weights are parame- 
terized by elliptic functions which satisfy the Yang-Baxter relations The model, 



hereafter referred to as the dilute Al model, has the outstanding property. For odd L's, 
the up-down symmetry of the local state is broken for nonzero elliptic nome . One may 
thus identify the magnetic field with elliptic nome. Actually, this identification leads to 
the desired critical exponent 5 = 15 of the Ising model for L = 3 [1^, |ll|. There are 



further supports for the conjecture that the dilute A3 model belongs to the same univer- 
sality class with the Ising model in a field, say, the central charge, critical exponents for 
bulk and a system with boundaries and so on p-[|l4 . 



From these two observations, one may naturally expect the Eg structure behind the 
dilute A3 model. As an evidence, we refer to the Thermodynamic Bethe Ansatz(TBA) 
approach in []T6|. These authors deal with the one dimensional counterpart to the dilute 
Al model, and discuss the finite temperature problem. They claim the dominating nine 
patterns (strings) of Bethe ansatz roots for L = 3 in the thermodynamics limit, and 
show that the string hypothesis leads to the Eg structure for the TBA equation. This 
observation originates from the numerical investigation near the zero elliptic nome. These 
strings are, however, quite involved and could lead to numerical difficulty with increase 



in the elliptic nome It may be preferable to have the derivation of Es structure 

free from specific choice of strings. 

In this report, we analyze the same problem through a different route based on the 
quantum transfer matrix (QTM) p3|, p^. The QTM approach has been developed recently 
as an alternative method in studies of ID quantum systems at finite temperatures 



The formulation yields the expression of the free energy by the largest eigenvalue 



of the object called QTM. The method has some advantages over the standard approach 
|T|-|^ based on the string hypothesis. For example, it yields a transparent derivation of 



the correlation lengths at finite temperatures ||4y, Q Once combined with the Yang- 
Baxter integrability structure, the evaluation of physical quantities at T > reduces to 
the study on analyticity of some auxiliary functions |31|,[^-|^. The crucial point in 



this approach lies in the choice of the auxiliary function, which is not necessary unique. 
One may adopt convenient functions having nice analytic properties. Here we choose 
them utilizing a family of commuting transfer matrices which includes QTM as the most 
fundamental one. We call them fusion QTMs. The functional relations among them allow 
the explicit evaluation of their eigenvalues. This strategy has been successfully applied 
to several models [^, There fusion QTMs come into play which corresponds to 

Yangian modules specified by only Young tableaux of rectangle shape. For the present 
model, this is no longer true. They constitute closed functional relations, however, lack 
nice analytic property. The absence of the property prevents the derivation of integral 
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equations (TBA). From trials and errors, we find it necessary to introduce variants of 
QTMs corresponding to "skew sliape" Young tableaux. A quantum analogue of Jacobi- 
Trudi formula [|6, J^, 5^, 50 1 as well as special features at "root of unity" play fundamental 
roles in the derivation of the functional relations. Numerical investigations on analyticity 
of "modified" QTMs show a remarkable coincidence of imaginary parts of locations of 
zeros with the exponents characterizing massive particles in scattering theory |51 



This validates a modest assumption on the existence of a strip in the complex plane where 
"modified" QTMs are free from zeros and singularities. Based on these, we successfully 
recover the TBA equation from the functional relations, namely, totally independent of 
string hypothesis. 

This paper is organized as follows. In the next section, we give a brief review on the 
dilute Ai models. A minimal guide to the QTM approach is sketched in section 3. The s/3 
fusion structure found in |52| will be explained in section 4. Based on these preparations, 
we introduce fusion QTMs parameterized by skew Young tableaux and their variants in 
section 5. The latter are found to satisfy closed functional relations. Section 6 is devoted 
to the observation of the analyticity of these functions. Assuming Conjecture 1 and 2, the 

TBA equation is derived from the functional relations. We conclude the paper with 
brief summary and discussion in section 7. 



2 dilute Ai model 

The dilute Al model is proposed in as an elliptic extension of the Izergin-Korepin 
model |5^. The model is of the restricted SOS type with local variables G {1, 2, ■ ■ ■ , L}. 
The variables {a,b} on neighboring sites should satisfy adjacency condition, |a — 6| < 1. 
The solvable weights are given by, 

a a ^ 9i{6 - u)9i{3 + u) 9i{u)9ii3 ~ u) 

a^a 9i{6)9i{3) 9,{6)9,{3) 

g4(2a-5) Sa.l9^{2a + 5) ^ 

V Sa ^4(2a + l) Sa ^4(2a-l)^' 

a±l^a ^ ar-^a _ 9i{3 - u)94{±2a + 1 - u) 

a a^a±l 9,{3)9^{±2a + 1) 

a „a a^a±l f Sa±i\m9i{u)9i{±2a - 2 + u) 
\u\ = \u\ - ' 



a, — I ( 
a ( 



a±l^a a^a ^ Sa ^ ^i(3)^4(±2a + 1 



a ± 1 a ± 1 ^ a ± 1 f0,{±2a + 3)9,{±2a - l).i/2 9^{u)9,{3 - u 



a±l a ^ a ^ ^|(±2a + l) ^ ^1 (2)^1 (3) 

a±l^ a 9^{2-u)9^{3-u) 



a ^a^l ^1(2)^1(3) 
a ^a^l _ (Sa-iSa+i^i/29i{u)9i{l-u) 



a±l^ a ^ SI ^ 9i{2)9i{3) 

a r^a±l 9^{3 - u)9^{±Aa + 2 + u) 
\u\ 



a±l^ a ^i(3)^i(±4a + 2) 
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Sa 6'i(3)6'i(±4a + 2) 
ei{3 + u)ei{±Aa - A + u) 



^i(3)^i(±4a-4) 
5aTi^i(4) 9ii±2a - 5)^9i{u)9i{±Aa - 1 + u) 



5A(2) ^4(±2a + l)^ 0i(3)0i(±4a-4) 
Here 6'i^4(x) = '(9i 4(Ax, r), 

oo 

^9i(x,r) = 2gi/Smxn(l-2g^"cos2x + g^")(l-g2^ 

n=l 

oo 

^4(x,r) = n(l-2g'""'cos2x + g^"-2)(l-g2"), 



otherwise .(11 



71=1 



and q = exp(— r). A is a parameter of the model specified below and Sa denotes 

Q - ( naMM 

The model exhibits four different physical regimes depending on parameters, 

• regime 1. < m < 3, A = 4(f^, L>2 

• regime 2. < m < 3, A = f[^^, L > 3 

• regime 3. 3-f<M<0, A = fjgfj, L > 3 

• regime 4. 3-f<M<0, A = j^j^y L>2 



We are interested in regimes 2 and 3. The central charges and the scaling dimensions of 
the leading perturbation have been evaluated in [1^, |lT|. Specializing to L = 3, they 
read 



regime 2. c = ^, A 



16 

regime 3. c = |, A = y|- 



Obviously, in connection with the Zamolodchiklov's argument, regime 2 attracts more 
attention. We will see, however, that both of them can be treated on a same footing. 

The one particle excitations are examined in |T5[. Eight particles are identified and 
their masses are summarized by a single formula. 



rrii = > sm 
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where 
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set of allowed a's 


1(1) 


{1,11} 


2(3) 


{2,10,12} 


3(5) 


{3,9,11,13} 


4(7) 


{4,8,10,12,14} 


5(8) 


{5,7,9,11,13,15} 


6 (6) 


{6,8,12,14} 


7(2) 


{7,13} 


8(4) 


{6,10,14} 



Table 1. 



A number k in the bracket means that it corresponds to the k— th light particle. These 
exponents appear in many contexts, diagonal scattering matrices, Weyl reflections with 
respect to Coxeter element and so on [^. They will re-appear in a novel context later. 



3 Quantum Transfer Matrix 

The row to row transfer matrix Tj^tj^(m) is defined by 

M 

(TRTRH)ffi = n^'0!^' 



+1 



From the commutativity of Tj^trI"^) with different spectral parameters, it is natural to 
define the Hamiltonian of associated ID quantum chain by 

d 

-^lnTRTR(M)k=o 



as in |jT6[. Here e = —1, (1) labels regimes 2 (3). In order to evaluate the free energy , 
we adopt the method of the quantum transfer matrix which is defined in the following 
"staggered" manner. 



(TQTMK^))S=n 

j=l a,2j-\ 





b2j 


0.2 j + l 




u+iv 






u — iv 




a2j 


a2j 





b2j+i 
b2j 



Note that the number of sites A^, sometimes referred to as the Trotter number, has no 
relation to the real system size M and is even by the construction. We further remark the 
commutative property of QTMs, 

[TQTM{u,v),Tq^M{u,v')] = 0. 

The free energy per site is represented only by the largest eigenvalue of Tqxm 
and u = 







'TV' 



lim —InTr exr>(—dHe) 

M-,oo M l-K y 



— lim ln(^the largest eigenvalue of 7qxm(^ = ~^iV' ^ 



0) . 
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Explicitly, the eigenvalue Ti{u,v) of TQ^y[{u,v) is given by 



2 2 (^(f — l/2z) 2 2 (^(i; — l/2z) + 1/2?) 

^ -U. 3 1 Qiv + 5/2t) 



2 '"^^ 2 'g(v + l/2i)' 

Ar/2 

g(t;) := Y[h[v-Vj] 

+ (3/2 - - (3/2 - u)i]^N/2 

h[2i]h[3i] 



h[v] ■.= 9^itv), 



where w = exp^ini/ (L + 1) (£ = 1 for the largest eigenvalue sector). The parameters, {vj} 
are solutions to "Bethe ansatz equation" (BAE), 

0(i;,+z) _ Q(v^-i)Q{v^ + 2i) N 

(P{vj-t) Qiv,+i)Qiv,-2i)' ' ' 2' ^ ' 

Now the difficulty arises. comes into the expression of the free energy through u, as well 
as the rhs of (0). Thus a simple-minded transformation of BAE into an integral equation 
using the root density function is not legitimate in contrast to T = problems. A naive 
extrapolation to A^ — > oo by numerics may bring about errors. One must thus devise 
other methods in encoding the information of BAE roots. In this report, we make use of 
the commuting family of the transfer matrices (fusion hierarchy). Indeed, it has been 
known that the functional relations among them, which hold for arbitrary A^, can be an 



alternative to BAE |54, In the next section, we shall argue the s/3 fusion structure in 
the dilute model. 



(3) 



4 sh fusion structure 



The s/3 type fusion structure in the dilute model has been discussed in details ^ . 



This comes from the singularity of the RSOS weights at -u = ±3 where the face oper- 
ator becomes projectors. A desired subspace can be picked up from tensor products of 
spaces by using these projectors. (For explicit procedure, see[0.) The adjacency matrices 



are identified with Young diagrams, and combinatorics of tableaux describe their tensor 
decomposition. 

We are interested in eigenvalues of fusion QTMs. In this view point, the most relevant 
fact is that these eigenvalues are again expressible in terms of "Young tableaux" depending 
on spectral parameters. Let three boxes with letters 1,2 and 3 represent the three terms 
in eigenvalue of the quantum transfer matrix; 

Ti(«,t;)=a + a+a- 



^ Results on other choices of auxihary functions for several models, see pn, pa, BGL 
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Obviously, the eigenvalue of a fusion QTM can be represented by sum over products of 
"boxes" with different letters and spectral parameters. The point is that the assembly of 
such boxes can be identified with semi-standard Young tableaux (SST) for sl^. We present 
a simple example. (See ^ ^ |50|, ^ for details.) By the fusion procedure, one 



can construct a transfer matrix of which auxiliary space acts on a symmetric subspace of 
\/x\/. The set of the SST 



«2 



{ii < {2) is associated with this subspace. The eigenvalue 



of the transfer matrix is then represented by, 

«1<«2 



V+i 



(5) 



Similar construction leads to fusion models based on general Young tableaux, of which 
eigenvalues are expressed by their shapes. On each diagrams, the spectral parameter 
changes +2i from left to right and —2i from top to the bottom. 

We restrict ourselves to rectangular shapes in this section. There are still interesting 
properties, some of which are crucial for the discussion in the next section. First, due to 
identities, 



V+l 

v—i 



v+i 
v—i 



5 5 
(f + -i)(l)iv - -i)\3 



v+l 



v+2i 



v-2i 



3 



+ (2 " 



,9 
'2 



the QTMs from 2 x m ( 3 x m) Young tableaux reduce to those from 1 x m (or just 
scalars). Second, the eigenvalues of 1 x m fusion QTMs have the "duality" in the sense 
of eq.(|^). For the explanation, we introduce renormalized 1 x m fusion QTMs Tm{v) by 



Tm{v) 



fmiv] 



E 



ii 


i2 




is 



where the semi-standard (SS) condition ii < ^2 ^ ■ ■ ■ ^ is imposed on the summation. 
(From now on we suppress the dependency on u which must be kept identical for all 
fusion QTMs.) The spectral parameters, v — i{m — 1) ■ ■ - v + i{m — 1), are assigned from 
left to right. A renormalization factor, which is the common factor of the expressions from 
tableaux of length m, is given by 



m— 1 

fm{^) ■■= n 



2m - 1 



j))(p{v - i 



,2m - 1 



Then the resultant T^'s are all degree 2N w.r.t. h[v + shift ], and have a periodicity due 
to Boltzmann weights; 

Tmiv + 16/52) = Tm{v). 
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From the s/3 structure, together with the above property, the following functional relations 
are valid. 



Tm{v - i)Tm{v + i) = gm{v)Tm{v) + Tra+l{v)Tm^i{v), m > 1 

gmiv) = 0(t; + i(m + 3/2))0(t;-2(m + 3/2)), 
T-iiv) := 

To{v) := Mv). (6) 

The periodicity, (piv + 16/5i) = leads to gm+ieiv) = gm{v) , {gs+m{v) = gmiv + 

8/5i)) m > and g5-m{v) = gmiy ± 8/5i), (0 < m < 5). As the adjacency matrices are 
vanishing, Tq{v) = Tt{v) = 0. Thus one deduces the duality relations, 

8i 

Tm{v) = T^-rn{v + — ) , m = 0, ■ ■ ■ , 5 (7) 
5 

and Tjn+^iv) = T^iv + ^),m > —1 for the solutions to eq.(|^). These relations can 
be in principle proved by using explicit fusion weights. We have at least checked the 
validity numerically. Thus, for the dilute case, only 5 "one- row" fusion QTMs need 
consideration. 

Finally we discuss the closed set of functional relations among rectangular types, 
motivated by success of several other cases. From T— system (|^), we are led to introduce 
the "3^—" functions as in [|^, |7 



y„(.):=^^±l)^f^.(l<,„<4). 

gm[v)Tm[v) 

which satisfy closed set of equation, 

y,.(v - r)y^> + = + y.... W)(i + y-t")) ,1 < „ < 4). (s) 

1 + ym[v) 

where y^iy) = 0. 

In p^ , |38| , TBA equations are re-derived from such subsets of fusion hierarchy. 



as y equations possess good analytic properties; both sides of equations are analytic and 
nonzero and have constant asymptotic behavior as |f | ^00 (ANZC) in certain strips. By 
taking the logarithmic derivative of both sides, and applying Cauchy's theorem, one can 
directly rewrite y equations by coupled integral equations which are nothing but TBA 
equations. 

The situation is different for the present model; numerically, we find "intrinsic" zeros 
and singularities of 3^m's in the strips = [—1, 1] and '^v = [—1.6, — 1]U[1, 1.6]. Therefore 
the above argument can not be applied to (||) 

Summarizing, the restriction to rectangle shapes leads to a closed set of equations, but 
with poor analytic properties. 

In the next section, we thus consider a wider class of skew Young tableaux which does 
lead to finitely closed functional relations having the desired property. 



8 



5 quantum Jacobi-Trudi formula and fusion hier- 
archy 



Let n and A be a pair of Young tableaux satisfying fji > Aj, Vi. We subtract a diagram A 
from fi. The resultant "narrower" one, consisted of (yUi — Ai, /i2 — A2, ■ ■ ■) boxes is called a 
skew Young diagram fi — X. (The usual Young diagram is the special case that A is empty, 
and we will omit A in the case hereafter.) In the theory of symmetric polynomials, the 



Figure 1: An example of a skew Young Tableau, (4,4)-(2). In the Ihs, the Young diagram 
(4,4) with the shaded part (2) is given, while we have the corresponding skew diagram in 
the rhs. 

Schur function deserves the one of the fundamental objects. Generally, a Schur function is 
defined in correspondence to a skew Young diagram. The Jacobi-Trudi formula asserts its 
representation by a determinant of a matrix of which elements are given by those related 
to "one-row" diagrams or " one-column" ones. An analogous formula holds for the present 
situation EBI, ESI SI CTl. 



Consider a set of semi-standard skew Young tableaux of the shape /i — A. As remarked 
in the previous section, spectral parameters are assigned to boxes in a tableau in such a 
manner that they change by +2i from left to right and —2i from top to the bottom. We 
fix the spectral parameter of the "top-left" box by 1; -|- — where fi[ denotes the 
depth of the tableaux. One regards each box in a tableau as an expression under the rule 



Figure 2: The spectral parameter v + z(/ii — is assigned to the hatched place 

(§) with the shift of the spectral parameter. Then the product over all constituting boxes 
yields an expression for a tableau. 

Theorem 1 Let T^/xiv) he the sum of the resultant expressions over the SST divided by 
a common factor, Hjii /f(j-Aj(^ + ^(/^'i— /^i+yUj + Aj — 2j + l)). Then the following equality 
holds. 

%/\{v) = det i<j,k<f,[{Tf,j-Xk-j+k{v + iifJ^'i - fii + fij + \k - j - k + 1))) (9) 
where T^^o := 0. 
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We regard this as a quantum analogue of the Jacobi-Trudi formula. The proof utilizes the 
decomposition rules for products of two tableaux with spectral parameters^, |50| , 
which is quite to parallel to those for Young tableaux. One must only pay attention that 
the allowed positions of a box is restricted by its spectral parameter. T^/x{v) may be 
naturally identified with the eigenvalue of QTM corresponding to fusion fi — X. The proof 
of this is beyond the scope of the present report. The important fact for our purpose 
is that T^/x{v) defined in such a manner is an analytic function of v due to BAE, and 
contains Ti{v) as a special case. The former assertion is not obvious from the original 
definition by the tableaux, but it is trivial from the quantum Jacobi-Trudi formula. 

In the same spirit, we introduce A^/x{v), which is analytic under BAE, from T^j\{y) 
by putting Tm>&{y) = in the latter, 

K^/x{v) := %/x{v)/.{T^>^{v) ^ 0}. 

The pole-free property of K^ixiy) is apparent from (|]). 

The following "dualities" are simple corollaries of eq.(|^) and (||). 

Lemma 1 

(10) 

(11) 
(12) 

(13) 

(14) 

(15) 



(16) 

To establish a closed set of functional equations we need further 

Lemma 2 . 

The following relations hold, 

fiiv -2i)f2{y + i) 

A M = ^^(t^ + i) A(4,2)/(i)(t^-7i) 

(6,4,4)7(3,3) 1 /2(t; + i) (t){v - t5/2)<f){v - t3/2)<l){v + t7/2)<f){v + t9/2) ^ ' 

13 7 7 13 

A(l0,7,7,4,4)/(6,6,3,3)(^^) = ' TrO<^(^ ~ oO0(^ + ^^^(^ + ^(6,4,3)7(3,2) (^^) , (19) 



A(4,4)/(2) 


[v) 


2 

= A(3,1){V+ -i), 




A(4,2)/(l) 


iv) 


3 

= A(4,3)/(2)(t^ + -0, 




A(4,4)/(3) 


[v) 


A / 7., 

= A(4,i)(t; + 




A(7,4,4)/(3,3) 


iv) 


= A(4,4,l)/(3)(t^+ 


), 


A(6,4,3)/(3,2) 


iv) 


A / 8 

= A(5,4,2)/(3,l)(t^ + ^ 




A(7,7,4,4)/(6,3,3) 


iv) 


= A(7,4,4,i)/(3,3)(f 




A(7,7,4,4,l)/(6,3,3) 


iv) 


= A(io,7,7,4,4)/(6,6,3,3) 


8. 

iv + -^ 
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11 5 9 15 

A(3,i)(t; - lli)A(4,3)/(2)(^^ - (20) 

(0) can be shown trivially from SS condition. To prove the rest, it is convenient to go 
back to the original definition of T. For example, we consider (0). 

^6,4,4)/(3,3) (^^) = A(6,4,4)/{3,3)(^^) + ro(t;)To(f + 2i)T8(t;) 

f5{v + i)f3{v-5i)f2{v + 6i) 
f2{v + i)h[v + 5z)/4(v - Al) 
h{v + i) T^{v-hi)T2{v-m) 



f2{v + i) (f){v - i5/2)(j){v - i3/2)(t){v + i7/2)(f){v + i9/2) ' 



We have applied in the first line. The second line follows from SS condition; 3 boxes in 
a column reduce to a scalar, so that the configurations in other two rows are completely 
independent. 



Figure 3: Decomposition of the skew diagram (6,4,4)/(3,3). 



In the last equation, the periodicity brings about a new feature, namely, further de- 
composition of the products of T, 



Tsi^v - 5i)T2{v - lOi) = A(4,2)/(i)(^^ - 7i) + 



f2{v - lOt)Mv - 5t) 



T,{v-7i). 



Now that Ts{v) = T:^{v) = Tq{v + |z) as discussed in the previous section, the proposition 
reduces to 



To f To f + 2t)To{v + -t) = I / FT 



To{v -7i + -i), 



which is proved by representing Tq by 0. The other two relations are proved similarly. 

We define the following fundamental objects (referred to as "modified fusion QTMs" 
in the introduction); 



Definition 1 

T«(t;) 



A(i)(t^) 

7 8 

A(4,l)(^^ + — «)/0(^^ + -O 
A(7,4,4)/(3,3)(t^)/{0(^^ + 3/2i)(P{v - 3/2l)} 
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T(^)(t;) := A(5,4,2)/(3.i)(^^) 



9 7 8 9 

A(7,7,4,4)/(6,3,3)(f + ^^^/{^(^ + + + ^^)} 



13 ., 

A(4,2)/(l)(t^ + 
A(2)(^^+ 



(21) 



In terms of h[v+ shift ], they are of degree, 2A^, 3A^, AN, 5N, QN, AN, 2N, 3N, respectively. 
Proposition 1 The following "Eg type T-system" holds among them, 



<^(^ + ,_)T(2)(^)+To(t;±-z), (22) 
To(^;)To(^;±i)+T«(t;)T(3)(^;), (23) 
To{v ± ^t)n{v ± ^t) + T('\v)T('\v), (24) 

To(t;)To(t; ± U^v ± ^z) + T('\v)T('\v) , (25) 
5 5 

To(^ ± ^^)n{v ± i)To(^; ± + T(^\v)T'-'\v)T('\m 

To{v ± U)To{v ± ^0 + r(^)(t;)r(^)(t;), 
5 5 





[v — 


I 

10 ^ ^ 


1 , 
10 ^ 


2^(2) 


[v — 


10 ' ^ 


1 

10^'* 


2.(3) 


^"^ — 


^,)T(3)(^ + 
10 ^ ^ 


1 

10^'' 


rjn(A) 


[v — 


10 ^ ^ 


1 , 
10 ^ 


rpi5) 


[v — 


10 ' ^ 


1 

10^'* 


rjn(6) 


[v — 


10 ^ ^ 


1 , 

10 ' 




X — 


l^)r(^)(x + 


1 , 
10 ^ 


2.(8) 


[v 


lz)T(«)(t; + 


1 

10^'* 



To{v)To{v±li)+T<^''\v) 
5 



(27) 
(28) 



(29) 



Here and in the following we use g{x ± y) to mean g{x + y)g{x — y) for some function g. 

Note that the difference of the arguments of Ihs is much smaller (|i) than the difference 
{2i) in the original 5/3 type fusion equation. This leads to the analytical property shown 
in the next section. 

Let us sketch the derivation of the above equations. We consider the decomposition of 
T*^"^'s and their dual partners. The simplest case, (|22|), follows from consideration on the 



decomposition of Ti{v — 4i)T4(f + i). From (|^) we have, 

T,{v - Ai)T,{v + i) = A(4,i)(t;) + T,{v)To{v - 3i). 

Duality (|^) allows us to rewrite T/^iT^) by of Ti(Tq). After the shift v ^ v + 0.7i together 
with periodicity, (H) is estabhshed. Eqs. (||), (|1),(|2|) (|3) and (H) can be proved 
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similarly. Next we consider (|26| ) by utilizing, 



A(7,7,4,4,l)/(6,3,3)('i^)A(io,7,7,4,4)/(6,6,3,3)('y + 5i) 

= To{x ± i)To{x ± 9i)To{x + ni)T5{x ± 4i)T5{x ± 6i)T5{x + Ui) 

+A(l0,7,7,4,4,l)/(6,6,3,3)('i^ + ^'i')^{7,7,4A)/(e,3,3)iv + 0- 

The eqs.(^) and ( |TI1) allow us to rewrite the rhs by T^^\T^^^ and T^^\ The second factor 
in the Ihs can be rewritten by T^^^ after transformations A(io,7,7,4,4)/(6,6,3,3) A(6,4,3)/(3,2) 



((|19f) in Lemma 2) — A(5 4 2)/(3,i) ((M) in Lemma 1). The first term can also be reduced 
to T(^) with one further step ((|16|) in Lemma 1). By shifting v ^ v — O.li and dividing 
both sides by the factor 0(f ± ± ± |i)0(f ± 2i), we obtain (^61). (|29|) can be 

proved with the aid of (^), (^) and (|18D. 

6 ^8 Thermodynamic Bethe Ansatz from Y-system 

Now we transform the functional relations into "gauge invariant" forms. Let us define Y— 
functions by combinations of T^""^ ; 

Definition 2 



y(2)(^)y(4)(^) ,W4W„. _ Ti^KjoE^ 



then functional relations follow from the T— system. To write it down neatly, we assign 
numbers to nodes in the Dynkin diagram of Eg as in the following figure. 



Figure 4: The Dynkin diagram for Eg. 
We denote a ~ 6 if a and b are adjacent in the Dynkin diagram. 
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Theorem 2 Functional relations among Y— functions exhibit the structure in the 
following form, 

+ = + a = l,---,8. 

This coincides with the Eg case of the universal Y— system in In order to reach our 
final goal, the TBA equation, further information is required on analytic structures of 
), 1 + Y^^-^v), a = 1, ■ ■ ■ , 8. We employ numerical calculations for some fixed values 
of q, N and (3 for this purpose. Remark that only the largest eigenvalue sector needs 
examinations. This is a drastic simplification from standard string hypothesis argument 
in which analyses on all excitation sectors are, in principle, necessary. Though we have 
performed numerics for relatively small values of N, intriguing patterns are already found 
for zeros ofT^^^v). First, these zeros are symmetric with respect to real axis, which assures 
that T^"^) (f)'s are real on the axis. Second, imaginary parts of coordinates of zeros have 
the remarkable coincidence with the exponents in Table 1. See Table 2 for the example 
of the case u = —0.08, q = 0.3, N = 12. We summarize them as a conjecture for arbitrary 
N. 

Conjecture 1 Zeros ofT^'^^ distribute along approximately on the lines, ~ ±0.1(aj ■ + 
1). The set {aj} agrees with {a} for the particle j in Table 1. 

In addition, asymptotic behaviors are specified as. 

Conjecture 2 . F^"-* 's have bounded asymptotic behaviors in the largest eigenvalue sector 
of QTM. Explicitly, their limiting values read, 

r«(±oo) = 2(l + x/2), Y^^\±oo) = il + V2)i5 + 3V2), 

r(3)(±oo) = 6(l + y2)(3 + 2v^), F(^)(±oo) =4(4 + 3v^)(5 + 3y2), 

r{5)(±oo) = 3(3 + 2^2)2(5 + 4\/2), r(^)(±oo) = 4(2 + y2)(4 + 3^2), 

y(7)(±oo) = 5 + AV2, r(^H±oo) = 4(4 + 3^2). 

Combining these with the T— system, we have 

Lemma 3 . Assume that above conjectures are valid. Then y(")(f) and 1 + Y^"'\v) 
are Analytic, NonZero and have Constant asymptotic behavior (ANZC) in strips '^v G 
[—1, 1], [—0"^, 0"^], respectively. 

y(^){v) = F(")(t;) for a ^ 1 and, 

Y^){v) = Y^^\v){k{v + t{l + u))k{v - z(l + u))y 

where e = 1(— 1) for u > 0{u < 0) , u = u/10 and 

_ , '^2ii7rv/4:,T') ^N/2 
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Here t' = 8r is introduced so as to respect the periodicity of Y— functions along the real 
axis. Remark that in deriving Lemma 3 we actually do not need such fine structures as 
given in Conjecture 1. Lemma 3 is robust in this sense. 

Thanks to the simple identity k,{v + i)K,{v — i) = ±1, the Ihs of the theorem 2 can 
be replaced by Y^^^{v — i)Y^'^^{v + i). Then we are in position to apply the Cauchy 
theorem directly and transform logarithmic derivatives of resultant algebraic functional 
relations to a set of coupled integral equations. Note that the resultant equations depend 
on Trotter number only in the renormalization factor of Y^^^ through the combination, 
roughly speaking, uN. The limit — >■ oo is thus performed straightforwardly. 

Let us introduce Fourier transformations by 



^ /•Wt/tt a ^ ^ . ^ 

F{x)= e-''"'F{v)dv, = — Ve™^"F(x, 

J-16r/7r ZTT ^ 



with 6 = 7r^/(16r),x„ = n6. We also denote by lnF(a;) the FT of logarithm of F{v). 

Theorem 3 By Lemma 3, one can prove a set of nonlinear integral equations among 
lnY^^\v) for arbitrary N. 

By taking the limit N ^ oo, the resultant equations read in the Fourier space as, 

\nY^^\x) = -e6a,iM^)+Ca,t{^Ml + Y'^'^){x) 

s(x) = — 

2 cosh X 

Ca,b{^) = S{X){2I -C'^^)a,b, (30) 

where (3 = 207r/5, and C^^ denotes the Cartan matrix for E^. 

The free energy is expressed via Y— functions with the aid of (p2D. 

-/3/ = -/5eo-/36i*s(0) + s*ln(l + rW)(0) e = 1, (31) 
= /3eo-,96o*Po + <?*ln(l + (F(*))-')(0) e = -1 (32) 
eo := A[ln(^i(7r/16)i9i(37r/8))]' 

g sinh 6x b i ) ^^'^^ ~'~ ^^^^ ^'^^ 

sinhl6x' sinhl6x ' 

Po(x) := ag(x):=5(x)((l-C(a;))-i)i_,, (33) 

16 / 

where A*B := f_-[l)^ A{v-v')B{v')dv' . The final expressions are given in different forms 
depending on the magnitudes of Y in the vicinity of the origin , (|y*^"''| ^ 1(^ 1) for 
e = l(-l)). 

The set of nonlinear integral equations is nothing but the TEA equations obtained in 
16| , and the free energy also coincides with their result. (Note that the normalization of 



our hamiltonian is 207r times the one in 

Therefore, we have reached the same results, however, completely independent of the 
string hypothesis. 
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7 Summary and Discussion 



In this report we have apphed the recent developed QTM method to the ID quantum chain 
related to the dilute As model. Modified QTMs described by skew Young tableaux have 
played a role, which has not been observed in former successful applications of the QTM 
method. The Eg structure of the TBA equation has been recovered without hypothesis on 
specific forms of dominant solutions. This gives further supports and consistency of the 
"Trinity" among minimal unitary CFT theory p = 3 perturbed by 0i^2, the Ising model 
in a field and the dilute ^3 model. 

There may be a direct connection to E^. In the context of the analytic Bethe ansatz, 
the dressed vacuum form for E^ case , consisted of 249 terms, is already obtained explicitly 



as shortly remarked in It would be interesting to show the direct reduction of 
these 249 terms to only 3 terms. 

The method employed here is in principle applicable to dilute A4, Aq model of which 
underlying symmetries are expected to be of E-j and Eq. Partial evidences are found very 
recently in a different context [^. We hope to report the study on them, as well as the 
extensive numerical studies on the Eg case in the near future. 
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Table 2 





exponent 


locations of zeros 




1 


±0.0160536139 ± 0.194948867i, ±0.0538522756 ± 0.192860391i, 

±0.124085873 ± 0.179123531?: 




11 


±0.0161627786 ± 1.19901092?, ±0.0544016497 ± 1.19859179?, 
±0.129454048 ± 1.19555064i 




2 


±0.0160900147 ± 0.295944203?, ±0.0540354627 ± 0.294252469? 
±0.125806251 ± 0.282759557? 




10 


±0.0161506232 ± 1.09758994i, ±0.0543363518 ± 1.09656338?, 

±0.128545339 ± 1.08913466i 




12 


±0.0161758534 ± 1.30002265?, ±0.0544655721 ± 1.30003507?, 
±0.130002895 ± 1.30016934? 




3 


±0.0160996137 ± 0.396199839?, ±0.0540833902 ± 0.394610881?, 

±0.126249978 ± 0.383723324i 




9 


±0.0161379988 ± 0.997143234i, ±0.0542727134 ± 0.995931521i, 
±0.12791469 ± 0.987294006? 




11 


±0.0161709223 ± 1.19860006?, ±0.0544390469 ± 1.19800087?, 
±0.129624985 ± 1.19357709? 




13 


±0.0161774898 ± 1.40026185?, ±0.0544730978 ± 1.40037631?, 
±0.130044745 ± 1.40126236i 




4 
8 


±0.0161030779 ± 0.496280551?*, ±0.054100446 ± 0.49472356?*, 
±0.126398274 ± 0.484020597?* 

±0.0161379988 ± 0.896993334?, ±0.0542433172 ± 0.895721491?, 

±0.127657169 ± 0.886719531i 




10 


±0.0161605145 ± 1.09815223i, ±0.0543870319 ± 1.09736514i, 
±0.129124614 ± 1.09163793? 




12 


±0.0161742694 ± 1.29883901?, ±0.0544557686 ± 1.29834129?, 
±0.12979684 ± 1.29464801? 




14 


±0.0161781274 ± 1.50032712?, ±0.0544760749 ± 1.50046909?, 
±0.130062127 ± 1.50155269? 



20 



T'(5) 





±u.uioiU4yuo / ± u.oyDoiuyoDz , ±u.uo4iuyzoi ± u.oy4/Dozo/z , 

±0.126463632 ± 0.584119659?:* 




i 


ztU.UiOizy M ± U. /yoy44Uoo«, ±U.U04zc)zzyo/ ± U. /yoDOzoioZ, 
±0.127566253 ± 0.786537192i 




Q 


±0.128903306 ± 1.09163793^ 




11 


±0.0161643894 ± 1.19839097i, ±0.054406539 ± 1.19770481i, 
±0.129330606 ± 1.19269022i 






itu.uioi ( oooyy ± i.oyoyu4z^, ±u.uo44oiziO( ± i.oyo4oo(y?, 
±0.129847801 ± 1.39493192i 




io 


4-0 ni Ri 78Qc;77 -I- 1 c;QQa7a79.i' -un nc;/i/i77Qn/iQ 4- 1 c;QQaiQS7.)' 
itu.uiDi / ooo / / zt i.oyyo / o / zz, ztu.uo44 / / oU4o zt i.oyyoiyo/z, 

±0.130078989 ± 1.59928832i 


± ^ 


u 

Q 
O 


±0.126821874 ± 0.684758808i 

ltU.UiDi0oZ04Z zt U.oy / yoZ4iOZ, zlzU.U04o0i04i4 zt U.oy / Uo4iooi, 

±0.128823242 ± 0.89083195^ 




1 2 


±0 01616^632'^ ± 1 2984^609? ±0 0544127321 ± 1 2977971? 
±0.129390339 ± 1.29296911i 






±0.130033087 ± 1.49763549i 


T(7) 
1 ^ ' 


7 
( 


4-n ni fil QQi^S97 4- n 7Q7/IQfi8Q0') -l-H ni^/l9fiQfi1 1^7 4- n 7QR/I/IQ1 i^Q-i 

ztu.uiDioyooz / zt u. / y / 4yooozz, zizu.uo4Zoooio / zt u./yD44oioyz, 

±0.128192944 ± 0.788942255i 




13 


±0 0161693285 ± 1 39904538? ±0 054432565 ± 1 39863761z 
±0.129668227 ± 1.39563892 




6 


±0.0161211658 ± 0.696762829i, ±0.0541902581 ± 0.695400475i, 
±0.127219855 ± 0.685875782i 




10 


±0.0161593273 ± 1.09824048i, ±0.0543817684 ± 1.0974928i, 
±0.129119115 ± 1.09207157 




14 


±0.0161762757 ± 1.4990832i ±0.0544660786 ± 1.49868902i, 
±0.129915909 ± 1.49573858 



The locations of zeros of fusion QTMs for N — 12,q — 0.3, u — —0.08. The corresponding 
exponents of Table 1 are given in the second column. 

* For these cases, we have a problem in the convergence by the simple Miiller algorithm. 
These values must be thus taken to be rather approximate. 
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